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Abstract. In this article, I give a pedagogical introduction and overview of percolation 
theory. Special emphasis will be put on the review of some of the most prominent of the 
algorithms that have been devised to study percolation numerically. At the central stage 
shall be the real-space renormalization group treatment of the percolation problem. As 
a rather novel application of this approach to percolation, I will review recent results 
using similar real-space renormalization ideas that have been applied to the quantum 
Hall transition. 

1 Introduction 

Imagine a large chessboard, such as occasionally found in a park. It is fall, and 
all the master players have fled the cold a long time ago. You are taking a walk 
and enjoy the beautiful sunny afternoon, all the colors of the Indian summer in 
the trees and in the falling leaves around you. Looking at the chessboard, you 
see that some squares are already full of leaves, while others are still empty [1]. 
The pattern of the squares which are covered by leaves seems rather random. As 
you try to cross the chessboard, you see that there is a way to get from one side 
of the board to the opposite side by walking on leaf-covered squares only. This 
is percolation — nearly. 

Before continuing and explaining in detail what percolation is about, let me 
outline the content of this paper. In Sect. 2, I will review some of the most 
prominent and interesting results on classical percolation. Percolation theory is 
at the heart of many phenomena in statistical physics that are also topics in this 
book. Beyond the exact solutions of percolation in = 1 and d — oo dimensions, 
further exact solutions in 2 < d < oo only rarely exist. Thus computational 
methods, using high-performance computers and algorithms are needed for fur- 
ther progress and in Sects. 2.2 - 2.4, 1 explain in detail some of these algorithms. 
Section 3 is devoted to the real-space renormalization group (RG) approach to 
percolation. This provides an independent and very suggestive method of an- 
alytically computing results for the percolation problem as well as a further 
numerical algorithm. 

While many applications of percolation theory are mainly concerned with 
problems of classical statistical physics, I will show in Sect. 4 that the percolation 
approach can give useful information also at the quantum scale. In particular, 
I will show that aspects of the quantum Hall (QH) effect can be understood by 
a suitably generalized renormalization procedure of bond percolation in d = 2. 
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This application allows the computation of critical exponents and conductance 

distributions at the QH transition and also opens the way for studies of scale- 
invariant, experimentally relevant macroscopic inhomogeneities. I summarize in 
Sect. 5. 

2 Percolation 

2.1 The Physics of Connectivity 

From the chessboard example given above, we realize that the percolation prob- 
lem deals with the spatial connectivity of occupied squares instead of simply 
counting whether the number of such squares has the majority of all squares. 
Then the obvious question to ask is: How many leaves are usually needed in 
order to allow passage across the board? Since leaves normally do not interact 
with each other, and friction-related forces can be assumed small compared to 
wind forces, we can model the situation by assuming that the leaves are ran- 
domly distributed on the board. Then we can define an occupation probability 
p as being the probability that a site is occupied (by at least one leaf). Thus our 
question can be rephrased in modern physics terminology as; Is there a threshold 
value pc at which there is a spanning cluster of occupied sites across an infinite 
lattice? 

The first time this question was asked and the term percolation used was in 
the year 1957 in publications of Broadbent and Hammersley [2-4]. Since then a 
multitude of research articles, reviews and books have appeared on this subject. 
Certainly among the most readable such publications is the 1995 book by Stauffer 
and Aharony [5] , where also most of the relevant research articles have been cited. 
Let me here briefly stimmarize some of the highlights that have been discovered 
in the nearly 50 years of research on percolation. 

The percolation problem in d = 1 can be solved exactly. Since the number of 
empty sites in a chain of length L is {l—p)L, there is always a finite probability 
for finding such an empty site in the infinite cluster at i — > oo and thus the 
percolation threshold is Pc = 1. Defining a correlation function g{r) ex. exp(— r/^) 
which measures the probability that a site at distance r from an occupied site at 
belongs to the same cluster, we easily find g{r) = p^ and thus ^ = — 1/ Inp 
{Pc — p)~^- Thus close to the percolation threshold, the correlation length ^ 
diverges with an exponent u = 1. 

In d = 2, the percolation problem provides perhaps the simplest example of 
a second-order phase transition. The order parameter of this transition is the 
probability P{p) that an arbitrary site in the infinite lattice is part of an infinite 
cluster [6], i.e.. 



and /3 is a critical exponent similar to the exponent u of the correlation length. 
The distribution of the sites in an infinite cluster at the percolation threshold 
can be described as a fractal [7], i.e., its average size M in boxes of length A'' 




(1) 
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increases as {M{N)) oc N^, where D is the fractal dimension of the cluster. As 
in any second-order phase transition, much insight can be gained by a finite-size 
scaling analysis [7,8]. In particular, the exponents introduced above are related 
according to the scaling relation /3 = {d—D)u [5]. Furthermore, it has been shown 
to an astonishing degree of accuracy, that the values of the exponents and the 
relations between them are independent of the type of lattice considered, i.e., 
square, triangular, honeycomb, etc., and also whether the percolation problem is 
defined for sites or bonds (sec Fig. 1). This independence is called universality. In 
the following, we will see that the universality does not apply for the percolation 
threshold Pc- Thus it is of importance to note that Pc for site percolation on the 
triangular lattice and bond percolation on the square lattice is known exactly: 
Pc = 1/2. Especially the bond percolation problem has received much attention 
also by mathematicians [9] . 




Fig. 1. Site (loft) and bond (right) percolation on a square lattice. In site percola- 
tion, the sites of a lattice are occupied randomly and percolation is defined via, say, 
nearest-neighbor sites. In bond percolation, the bonds connecting the sites are used 
for percolation. The thin outlines define the 3 clusters in each panel. The solid outline 
indicates the percolating cluster 

For higher dimensions, much of this picture remains unchanged, although 
the values of Pc and the critical exponents change. The upper critical dimension 
corresponds to d = 6 such that mean field theory is valid for d > 6 with exponents 
as given in Table 1. 

Applications of percolation theory are numerous [10]. It is intimately con- 
nected to the theory of phase transitions as discussed in [8,11]. The connectivity 
problem is also relevant for diffusion in disordered media [7,12] and networks [13]. 
A simple model of forest fires is based on percolation ideas [14] and even models 
of stock market fluctuations [15] have been devised using ideas of percolation 
[16] . Percolation is therefore a well-established field of statistical physics and it 
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Table 1. Critical exponents v and /3 and fractal dimension D for different spatial 
dimensions. For a more complete list see [5]. 
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d = 6-e 
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4/3 
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0.68 


0.57 


1/2 + 5£/84 
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3.54 
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12/5 


2.8 




D{p > Pc) 


2 


3 


4 


5 





continues its vital progress with more than 230 publications in the cond-mat 
archives [17] alone. 

2.2 The Coloring Algorithm 

As stated above, there are only few exact results available in percolation in d > 2. 
Thus in order to proceed further, one has to use computational methods. 

The standard numerical algorithm of percolation theory is due to Hoshen 
and Kopelman [18]. Its advantage is due to the fact that it allows to analyze 
which site belongs to which cluster without having to store the complete lattice. 
Furthermore, this is being done in one sweep across the lattice, thus reducing 
computer time. 

At the heart of the Hoshen-Kopelman algorithm is a bookkeeping mechanism 
[5]. Look at the site percolation cluster in Fig. 1. Going from left to right and 
top to bottom through the cluster, we give each site a label (or color) as shown 
in Fig. 2. If its top or left neighbor is already occupied, then the new site belongs 
to the same cluster and gets the same label. Otherwise we choose a new label. In 
this way we can proceed through the cluster until we reach a problem in line 3 
as shown in the left column of Fig. 2. According to our above rule the new site, 
indicated by the bold question mark, can be either 2 or 4. This indicates that 
all sites previously labelled by 2 and 4 actually belong to the same cluster. Thus 
we now introduce an index Id(-) for each cluster label and define it such that 
the index of a superfluous label, say 4. points to the right label, viz. Id(4) = 2. 
Proceeding with our analysis into the 4th row of the lattice, we see in the center 
column of Fig. 2 that we again have to adjust our index at the position indicated 
in bold. Instead of labeling this site as 4, we choose 2 = Id(4). And consequently, 
Id(3) = 2. In this way, we can easily check whether a cluster percolates from top 
to bottom of the lattice by simply checking whether the label of any occupied 
site in the bottom row of the lattice has an index equal to any of the labels 
in the top row. Furthermore, in addition to the top row, we only need to store 
the row presently under consideration and its predecessor. Thus the storage 
requirement is linear in lattice size L and not L'^ as it would be if we were to 
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store the full lattice. Last, the algorithm can also give information about all 

clusters, say, if needed for a fractal analysis of the non-percolating clusters. A 
Java implementation of such a coloring algorithm can be found in Ref. [19]. 
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Fig. 2. Schematic description of the Hoshen-Kopelman coloring algorithm of the site 
percolation problem on a square lattice as shown in Fig. 1. • denotes an occupied site, 
the numbers denote the cluster labels. The horizontal lines bracket the current and the 
previous row 



2.3 The Growth Algorithm 



When we want to study primarily the geometrical properties of percolation clus- 
ters, another algorithm is more suitable which allows to generate the desired 
cluster structure directly. This algorithm is due to Leath [20] and works by a 
cluster-growth strategy. The idea of the algorithm is that we put an occupied 
site in the center of an otherwise empty lattice. Then we identify its nearest- 
neighbor sites as shown in Fig. 3. Next we occupy these sites according to the 
desired percolation probability p. We identify the new, yet undecided nearest- 
neighbor sites, occupy these again with probability p and repeat the procedure. 
The cluster continues to grow until either all sites at the boundary are unoc- 
cupied or the cluster has reached the boundary of the lattice. For p < pc, the 
growth usually stops after a few iterations, while for p > Pc, percolating clusters 
are generated almost always. Thus besides giving information about the fractal 
structure of the percolating clusters, the Leath algorithm can also be used to 
estimate the value of pc- A Java implementation of the Leath algorithm can be 
found in Refs. [21,22]. 



6 Rudolf A. Romer 



• ■ 
• 

• • 



o 


■ 


o • 


• 


• • 


• 


• • 
















o 


• 





o 


• 


• 


o 


• 


• 







• 


o 







Fig. 3. Schematic description of tiie first tiiree steps in tiie Leatii growtli algoritiim 
of the site percolation cluster of Fig. 1. •, o, and denote an occupied, empty, and 
undecided site 



2.4 The Frontier Algorithm 

As we have seen in the last section, the outer frontier of the percolation cluster 
is well defined. The fractal properties of this hull can be measured [23,24] and 
shown to yield D\i = 1.74 ± 0.02. This suggests yet another algorithm for the 
determination ofpc [25,26]: Generate a lattice with a constant gradient Vp of the 
occupation probability p as shown in Fig. 4. By an algorithm similar to the one 




Fig. 4. An example for gradient percolation. All sites at p = 1 {p = 0) arc occupied 
(empty). The 15 large circles correspond to occupied sites. The 8 light shaded sites 
belong to the frontier, whereas the 7 dark shaded sites are part of the interior of the 
cluster or belong to other clusters. The 5 small circles correspond to sites in the empty 
frontier. The dashed line indicates the frontier generating walk. Note that only the 14 
sites with solid circles have actually been visited. The 6 other circles arc shown here 
just for clarity and need not be generated. According to (2), we have Pc = (8 • 0.4375 + 
5 • 0.75)/13 = 0.55769 



used in computing the hull of the percolating cluster, one traverses the frontier 
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Table 2. Various current estimates of pc for site and bond percolation on different 
lattices in d = 2 [28-33]. Exactly known values are emphasized. The lattices are classi- 
fied according to their number of first, second and more nearest-neighbors. The upper 
index gives the corresponding number in the dual lattice 



lattice 




site Pc 


bond Pc 


3, 12^ 




0.807904 




4,6, 12 




0.747806 




4,8^ 




0.729 724 




6^ 


honeycomb 


0.697043 


0.652 703 


3, 6, 3, 6 


Kagome 


0.652 703 


0.524 4 


4* 


dice 


0.584 8 


0.475 4 


3,4,6,4 




0.621819 




4^ 


square 


0.592 746 


0.500 000 


3^6 




0.579 498 




3', 4, 3, 4 




0.550 806 




3^42 




0.550 213 




3« 


triangular 


0.500 000 


0.347296 



of the occupied cluster and determines which sites belong to it and which belong 
to the empty cluster [24 26]. Then a very reliable estimate for the percolation 
threshold can be computed as 

where iVo,e denotes the number of sites in the occupied (empty) frontier and Pco,ce 
is the mean height of the associated frontiers, respectively. Note that instead of 
actually generating the percolation lattice, the algorithm instead proceeds by just 
generating the sites needed for the construction of the frontier. Thus instead of 
dealing with, say, O(i^) sites as in the Hoshen-Kopelman and Leath algorithms 
for a lattice in d = 2, one only needs 0{L) sites in the present algorithm [27]. 
This reduces computer time and estimates for p^ in a large variety of lattices 
can be obtained with high accuracy [28-33] as shown in Table 2. 

3 Real-Space Renormalization 
3.1 Making Use of Self-Similarity 

As mentioned in Sect. 2, the transition at Pc corresponds to a second-order phase 
transition and the correlation length ^ is infinite. There is no particular length 
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scale in the system and all clusters are statistically similar to each other. This 
self-similarity [7] is at the bottom of the rcnormaUzation description just as for 
the fractal analysis of the percolation clusters. 

We may use the self-similarity in the following way. Let us replace a suitable 
collection of sites by super-sites and then study percolation of the super-lattice 
[34-38]. In general, the occupation probability p' of the super-lattice will be 
different from the original p. Furthermore, if the extent of the collection of sites 
in the lattice was b, then the super-lattice will have a lattice constant b. Thus 
^ = b£,' and with ^ oc \p — Pc\~'' , we find 



and consequently. 



b\p'-p,\-^ = \p-p. 



log 6 



(3) 



(4) 



As an example, let us consider the bond percolation problem on a square lattice 
[35 37]. Here wc replace 5 bonds by a super-bond in the horizontal direction as 
shown in Fig. 5. Summing all probabilities for a connected, horizontal super-bond 



p'(l-p) 



4p (1-p) 



2p'(l-p)' 



2p'(l-p)' 



4p'(l-p)' 



2p'(l-p/ 



Fig. 5. The possible combinations of bonds (thin lines) that lead to a super-bond (thick 
line) together with their respective probabilities for bond percolation on a square lattice 



as shown in Fig. 5, we find that 

p' =p^ + 5/(1 -p) + 8/(1 - pf + 2/(1 - pf. 



(5) 



At the transition, we have p' = p and thus (5) has the solutions p = 0, 0.5, and 1. 
The first and last solution correspond to a completely empty or occupied lattice 
and are trivial. The second solution reproduces the exact result of Table 2. From 
4), we compute v = 1.4274 which is already within 8% of the exact result 4/3. 
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Thus the real-space RG scheme gives very good approximations to the known 

results. But beware, it may not always be that simple: the reader is encouraged 
to devise a similar RG scheme for site percolation on a square lattice. 

3.2 Monte-Carlo RG 

The scheme of the last section is approximate since it cannot correctly handle 
situations like the one in Fig. 6. In order to improve, we can construct an RG 
scheme that uses a larger collection of bonds. The total number of (connected 
and unconnected) configurations in such a collection of n bonds is 2", putting 
severe bounds on the practicability of the approach for analytic calculations. 
However, the task is ideally suited for computers. On the CD accompanying this 
book, I include a set of Mathematica routines that compute the real-space RG 
for a d = 2 triangular lattice. 



Fig. 6. Although the original bonds (thin linos) arc not connected, the RG procedure 
outlined in the text nevertheless leads to two connected horizontal super-bonds (thick 
Unes) 



In 1980, von Klitzing et al. [39] found that the Hah resistance Ru of MOSFETs 
at strong magnetic field B exhibits a step-like behavior which is accompanied 
with a simultaneously vanishing longitudinal resistance R. This is in contrast to 
the classical Hall effect which gives a linear dependence of i?H on B. Even more 
surprising, the values of i?H at the transitions are given by universal constants, 
i.e., i^, where i is an integer. 

Since its discovery this so-called integer quantum Hall effect (IQHE) was 
studied extensively [40,41]. Besides semi-phcnomcnological models simply as- 
suming a localization-delocalization transition more general theories considered, 
e.g., gauge invariance [42], topological quantization [43], scattering [44] and field 
theoretical approaches [45]. 

4.1 Basics of the IQHE 

A simple understanding of the IQHE can be gained by considering the Hamilto- 
nian of a single electron in a magnetic field. 



4 The Quantum-Hall Effect 
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where A denotes the vector potential and the Hamiltonian has been rewritten in 

guiding center coordinates X = a;— C, y = and relative coordinates C, rj [46]. 
Here, Wc = ^ is the frequency of the classical cyclotron motion and = {-^y^"^ 
is the radius of the cyclotron motion. The spectrum of this Hamiltonian is simply 
the harmonic oscillator with £"„ = (r?, + i) hujc, n = 0, 1, . . .. These Landau 
levels are infinitely degenerate since the Hamiltonian no longer contains X and Y. 
Thus the spectrum consists of 5-function peaks as indicated in Fig. 7. Introducing 
disorder into the model by adding a smooth random potential V{r) in (6) results 
in drift motion of the guiding center 

perpendicular to the gradient of V{r) (see Fig. 8). Furthermore, the degeneracy 
of the Landau levels is lifted, the ^-function density of states broadens [40], giving 
rise to a band-like structure as shown in Fig. 7. If the sample is penetrated by a 
strong magnetic field, the cyclotron motion is much smaller than the potential 
fluctuations. Consequently, the electron motion can be separated into cyclotron 
motion and motion of the guiding center along equipotential lines of the energy 
landscape. 




Fig. 7. Density of states (DOS), transversal and longitudinal conductivity as a function 
of Ef or, equivalently, filling factor z/f or [41]. The peak in the middle of the band 
represents one (5-function peak of the clean Landau model. Dark shaded regions of the 
density of states correspond to localized states. The thin dashed line (with non-zero 
slope) for axy indicates the classical Hall result 



The IQHE can then be understood as follows: assume that the center of the 
broadened Landau levels contain extended states that can support transport, 
whereas the other states are spatially localized and cannot. This is similar to 
the standard picture in the theory of Anderson localization [7,12]. Changing the 
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Fermi energy Ep or the filling factor Uf = 2Trl%pe = 2TrUpe/eB oc Ep, where pe 
denotes the electron density, we first have Ep in the region of localized states 
and both axx and axy are 0. When Ep reaches the region of extended states, 
there is transport, a^x is finite and axy = e^/h. Next, Ep again reaches a region 
of localized states and axx drops back to until we reach the extended states in 
the next Landau level. 

This picture suggests the following effective classical high-field model [47] of 
the IQHE: Neglecting the cyclotron motion (i.e., large B) and quantum effects 
(i.e., only one extended state) the classical electron transport with energy 
through the sample only depends on the "height" of the saddle points in the 
potential energy landscape V{r). One obtains a classical bond-percolation prob- 
lem [5], in which saddle points are mapped onto bonds. A bond is connecting 
only when the potential of the corresponding saddle point equals the energy of 
the electron Ep. From percolation theory follows [5] that an infinite system is 
conducting only when Ep = {V). Using this model one could already describe 
the localization-delocalization transition and thus the quantized plateaus in re- 
sistivity observed in IQHE [40] . But for bond percolation the correlation length 
diverges at the transition with an exponent of v = 4/3 which is in contrast to 
the value found in the QH experiments. 



4.2 RG for the Chalker-Coddington Network Model 

The Chalker-Coddington (CC) network model improved the high-field model 
by introducing quantum corrections [48], namely tunneling and interference. 
Tunneling occurs, in a semiclassical view, when electron orbits come close enough 
to each other and the electron cyclotron motions overlap. This happens at the 
saddle points, which now act as quantum scatterers connecting two incoming 
with two outgoing channels by a scattering matrix as shown in Fig. 8. Similar 
to bond percolation a network can be constructed such that the saddle points 
are mapped onto bonds. While moving along an equipotential line an electron 
accumulates a random phase which reflects the disorder of V{r). Results for 
this quantum percolation also show one extended state in the middle of the 
Landau band. The critical properties at the transition, especially the value of 
the exponent v ~ 2.4 [49], agree with experiments [50,51]. 

As explained for the bond percolation problem we now apply the RG method 
to the CC model. The RG structure which builds the new super-saddle points is 
displayed in Fig. 9. It consists of 5 saddle points drawn as bonds. The links (and 
phase factors) connecting the saddle points are indicated by arrows pointing in 
the direction of the electron motion due to the magnetic field. Each saddle point 
acts as a scatterer connecting the 2 incoming /i,2 with the 2 outgoing channels 

Ol,2 




with reflection coefficients and transmission coefficients ti, which are assumed 
to be real numbers. The complex phase factors enter later via the links between 
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Fig. 8. Left: Schematic plot of a smooth random potential V{r) with equipotential 
lines at _B = (V) indicated in black. Right: Equipotential lines of the same potential 
for E = (V) — i5max/2, (V^), and (V) +-Emax/2 corresponding to long dashed, solid and 
short dashed lines. Note the solid line percolating the system from top to bottom as 
indicated by the arrows 



Fig. 9. Left: A single saddle point (circle) connected to incoming and outgoing cur- 
rents li, Oi via transmission and reflection amplitudes t and r. Right: A network of 5 
saddle points can be renormalized into a single super-saddle point by an RG approach 
very similar to the bond percolation problem of Sect. 3. The phases are schematically 
denoted by the (jj's 

the saddle points. By this definition - including the minus sign - the unitarity 
constraint tf + rf = 1 is fulfilled a priori. The amplitude of transmission of the 
incoming electron to another equipotential line and the amplitude of reflection 
and thus staying on the same equipotential line add up to unity - electrons do 
not get lost. 

In order to obtain the scattering equation of the super-saddle point we now 
need to connect the 5 scattering equations according to Fig. 9. For each link the 
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amplitude of the incoming channels is defined by the amplitude of the outgoing 
channel of the previous saddle point multiplied by the corresponding complex 
phase factor e'*'=. This results in a system of 5 matrix equations, which has to 
be solved. One obtains an RG equation for the transmission coefficient t' of the 
super-saddle point [52] analogously to Eq. (5), 

depending on the products fi.j = U - . . .-tj, ri,,,j ^ Vi - . . .-Vj of transmission and 
reflection coefficients ti and of the i = 1, . . . , 5 saddle points and the 4 random 
phases (pk accumulated along equipotentials in the original lattice. For further 
algebraic simplification one can apply a useful transformation of the amplitudes 
U = (e^* + l)~^/^ and = (e~^'-|-l)~"^/^ to heights Zi relative to heights Vi of the 
saddle points. The conductance G is connected to the transmission coefficient t 
by G= |i|2eV/i [53]. 



4.3 Conductance Distributions at the QH Transition 

For the numerical determination of the conductance distribution, we first choose 
an initial probability distribution Pq of transmission coefficients t. The distribu- 
tion is discretized in at least 1000 bins. Thus the bin width is typically 0.001e/\//i 
for the interval t G [0, e/Vh]- 

Using the initial distribution Po{t), we now randomly select many different 
transmission coefficients and inscirt them into the RG equation (9). Further- 
more, the phases j ~ f, . . .4 are also chosen randomly, but according to a 
uniform distribution e [0, 2tt]. By this method at least lO'' super-transmission 
coefficients t' are calculated and their distribution Pi{t') is stored. Next, Pi is 
averaged using a Savitzky-Golay smoothing filter [54] in order to decrease sta- 
tistical fluctuations. This process is then repeated using Pi as the new initial 
distribution. 

The iteration process is stopped when the distribution Pj is no longer distin- 
guishable from its predecessor Pi-i and we have reached the desired fixed-point 
(FP) distribution Pc{t)- However, due to numerical instabilities, small deviations 
from symmetry add up such that typically after 15-20 iterations the distributions 
become unstable and converge towards the classical FPs of no transmission or 
complete transmission similar to the classical percolation case. Figure 10 shows 
this behavior for one of the RG iterations. The FP distribution Pc(G) shows a 
fiat minimum around G = O.be^/h and sharp peaks at G = and G = e^/h. 
It is symmetric with (G) = 0A98e^/h. This is in agreement with previous the- 
oretical [56,57] and experimental [58] results whereas our results contain much 
less statistical fluctuations. Furthermore we determine moments ((G— (G))'") of 
the FP distribution Pc{G). As shown in Fig. 10 for small moments up to m = 6 
our results agree with the work of Wang et al. [55], who computed moments 
m < 8.5. But more interesting is the fact that the obtained moments of the FP 
distribution can hardly be distinguished from the moments of a simple constant 
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Fig. 10. Conductance distribution at a QH plateau-to-plateau transition. The squares 
correspond to the fixed-point distribution, dashed and dot-dashed lines to the initial 
distribution and an unstable distribution, respectively. The solid line indicates a fit of 
the FP distribution Pc{t) by three Gaussians. Inset: Moments of the FP distribution 
Pc{G). The dashed lines indicate various predictions based on extrapolations of results 
for small m [55]. The dotted line denotes the moments of a constant distribution 

distribution thus indicating the influence of the broad flat minimum of the FP 
distribution around G = 0.5e^//i. 

For the determination of the critical exponent, we next perturb the FP dis- 
tribution sHghtly, i.e., we construct a distribution with shifted average Gq. Then 
we perform an RG iteration and compute the new average Gi of Pi{G). Tracing 
the shift of the perturbed average (?„ for several initial shifts Go, we expect to 
flnd a linear dependence of G„ on Go for each iteration step n. The critical ex- 
ponent is then related to the slope dG„/dGo [59]. Figure 11 shows the resulting 
V in dependence on the iteration step and thus system size. The curve converges 
close to u K. 2.4, i.e. the value obtained by Lee et al. [49]. Note that the "sys- 
tem size" is more properly called a system magnification, since we start the RG 
iteration with an FP distribution valid for an infinite system and then magnify 
the system in the course of the iteration by a factor 2". 

5 Summciry and Conclusions 

The percolation model represents the perhaps simplest example of a system 
exhibiting complex behavior although its constituents - the sites and bonds - 
are chosen completely uncorrected. Of course, the complexity enters through the 
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Fig. 11. Critical exponent z/ as a function of magnification factor 2" for RG step n. 

The dashed line shows the expected result i' = 2.39. Inset: The shift of the average Gn 
of P{G) is linear in Go- The dashed lines indicate linear fits to the data 



connectivity requirement for percolating clusters. I have reviewed several numer- 
ical algorithms for quantitatively measuring various aspects of the percolation 
problem. The specific choice reflects purely my personal preferences and I am 
happy to note that other algorithms such as breadth- and depth-first algorithms 
[60] have been introduced by P. Grassberger in his contribution [61]. 

The real-space RG provides an instructive use of the underlying self-similarity 

of the percolation model at the transition. Furthermore, it can be used to study 
very large effective system sizes. This is needed in many applications. As an 
example, I briefly reviewed and studied the QH transition and computed con- 
ductance distributions, moments and the critical exponent. These results can be 
compared to experimental measurements and shown to be in quite good agree- 
ment. 
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